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We study the electronic properties of h-BN/graphene/h-BN ABC-stacked trilayer systems using 
tight binding and DFT methods. We comment on the recent work of Ramasubramaniam et al. 
(arxiv:1011.2489) whose results seem to be in disagreement with our recent calculations. Detailed 
analysis reaffirms our previous conclusions. 
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Following our recent paper [lj, Ramasubramaniam et 
al. have investigated the electronic properties of h- 
BN/graphene/h-BN sandwiched structure and general- 
ized it on the configurations containing bilayer graphene 
(BLG) 0. They found, however, that their DFT calcu- 
lations are not consistent with our results [l[ obtained 
by means of combined tight binding (TB) and ab initio 
methods. While overall valuable, we believe that this ar- 
ticle does not invalidate our model and conclusions. The 
purpose of this paper is to explain the contradictions and 
to clarify different interpretations. 

The previous analysis of the electronic properties of 
graphene/h-BN structures [l|, |3[ was based on the mini- 
mal tight binding model with parameters determined by 
density functional theory (DFT) calculations. As far 
as we know, the basis containing localized orbitals is 
quite reasonable in calculations for materials consisting 
of boron, carbon and nitrogen. Tight binding method 
allows to obtain correct spectra of graphene and h-BN 
sheets in the whole Brillouin zone when three nearest 
neighbors are included in the model, while the qualita- 
tive description of dispersion relations in the vicinity of 
K/K' points is possible even when only nearest neighbor 
model is considered. We found the small discrepancies in 
the TB curve fitted to DFT results [l[ , but this fact can- 
not significantly influence the properties of the system, 
especially because the Dirac cones are preserved in DFT 
band structure. 

It has been assumed in the TB model that the positions 
of atoms within each layer are fixed: the same geometry 
was used in DFT calculations, i.e. no relaxation was al- 
lowed. For completeness, we have calculated the spectra 
of relaxed structure using exactly the same parameters 
as in the work of Ramasubramaniam (see also Methods 
in Ref. [2|). It is clear that the small differences in lattice 
structures do not strongly affect the n bands, since the 
Dirac cones are still preserved in the ABC trilayer con- 
taining N-C dimers. In addition, we would like to stress 
that DFT calculations described in Refs 0, Q were per- 
formed under almost the same assumptions as those for 
graphene on h-BN substrate [4j and the results given in 
Ref. [1] have been well reproduced. Therefore it can be ex- 
pected that the same parameters will give correct results 



for bilayer and trilayer graphene/h-BN configurations. 

The next important issue is to explain the statement 
in the footnote 26 in Ref.Qj. Figure 3 in Ref. [l[ shows 
the variation of the band gap with the effective value of 
energy difference between the layers. The effective value 
means exactly the U parameter in the TB hamiltonian 
(equation (1) in Ref. [1|) and it is not equivalent to the 
field in the inputs of DFT calculations, since the screen- 
ing effects are not included in TB model. The value of U 
related to each DFT band structure in Refs [H, Q is esti- 
mated by analysis of band shifts which are predicted by 
TB model: the shift of bands in the spectra determines 
the exact value of parameter U used in TB hamiltonian. 
This means that the value 5 V/nm in Ref. Q is equiva- 
lent to about 1.15 eV in the scale of Fig. 3 of Ref.[l[. It 
can be easily found that for this value of U the predicted 
band gap is about 40 meV , which is in agreement with 
the results of Ramasubramaniam shown in Fig. 5 (e) of 
Ref. @. 

Moreover, it is explicitly stated in the last section of 
Ref. [l[ (illustrated in Fig. 3 therein, where the system 
is compared to the freestanding BLG), that, in practice, 
the opening of theoretically predicted band gap of about 
230 meV may be very difficult to achieve, because, for 
this system, the very large or even unphysical electric 
fields are needed. It is also worthwhile to note that simi- 
lar highly symmetric ABC configuration containing B-C 
dimers was rejected as less promising, because the band 
gap can be at most 50 meV and cannot be larger even 
theoretically. Therefore, the band gaps in both B-C and 
N-C ABC trilayers hardly depend on the applied electric 
field. On the other hand, the most important feature of 
the considered trilayers is the presence of the Dirac cones 
in the vicinity of the Briluoin zone corners and the fact 
that the band gap can be created. 

Finally, contrary to the statement of authors of Ref. , 
we observe that, in the band structures shown in Figs 5 
(a)-(c) (ABA configurations) and Fig. 5 (f) (ABC with 
different dimers) in Ref.Q, the linear dispersion is not 
preserved near K/K' points due to the symmetry break- 
ing between graphene sublattices. It is typical of the 
graphene/h-BN systems that the mass of Dirac fermions 
increases with the value of a band gap [5(; thus in all 
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FIG. 1: Band structure of the relaxed ABC h- 
BN/graphene/h-BN trilayer containing N-C dimers. 

cases where the band gap is present the 7r bands should 
be parabolic. The effective mass is finite, which is a fun- 
damental difference when compared to the freestanding 
monolayer graphene. We have performed the calculations 
for the case of relaxed structures shown in Figs 5 (d)-(e) 
of Ref . Q and we found no gap between valence and con- 
duction bands (Fig. 1). Therefore, the statement that 
the band gaps are in range of 10-100 meV is not clear to 
us. 

It should be stressed that, indeed, the bilayer graphene 
inserted between the h-BN sheets can provide a system 
of immediate practical relevance for graphene - based 



devices, especially in view of recent experiments [6j. 
The mechanism of band gap tuning in h-BN/BLG/h-BN 
sandwich is very similar to that in the freestanding BLG 
and the sensitivity to the stacking is not as important 
as in the configurations containing monolayer graphene 
(graphene is randomly stacked on bulk h-BN 6]). On 
the other hand, the band structure manipulations in the 
ABC systems with single layer graphene are unique and 
could offer new possibilities of unusual applications. Un- 
fortunately, the exact ABC stacking is needed for the full 
realization of band gap tuning in this system. The lat- 
tice constants should be precisely matched, which is not 
certain until now pj. 
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